


f 


j! S nN 
iby 4 ‘ ts Me , 









Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1996-12 


Experimental performance studies of a plate 
heat exchanger 


Plath, Darren R. 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/8375 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 





LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


http://www.nps.edu/library Monterey, California USA 93943 


NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 





THESIS 


EXPERIMENTAL PERFORMANCE 
STUDIES OF A PLATE 
HEAT EXCHANGER 
by 
Darren Plath 


December, 1996 


Thesis Advisor: Ashok Gopinath 





Approved for public release; distribution is unlimited. 


Thesis 
P6123 





DUDLEY KNOA LIS 
NAVAL POSTGRADUATE SCHOO! 
MONTEREY CA 93943-5101 





REPORT DOCUMENTATION PAGE /—ramamoesonere our | 


Public reporting burden for this collection of information is estimated to average | hour per response, including the time for reviewing instruction, searching existing data 
sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any 
ollier aspect of this collection of information, including suggestions for reducing this burden, ty Washington Headquarters Services, Divectorace for Information 






Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction 
Project (0704-0188) Washington DC 20503. 


1. AGENCY USE ONLY (Leave blank) ah REPORT DATE 3. REPORT TYPE AND DATES COVERED 
December 96 Master’s Thesis 


4. TITLE AND SUBTITLE EXPERIMENTAL PERFORMANCE STUDIES . FUNDING NUMBERS 
OF A PLATE HEAT EXCHANGER 


6. AUTHOR(S) Darren Plath 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 


8. PERFORMING 
Naval Postgraduate School ORGANIZATION 
Monterey CA 93943-5000 REPORT NUMBER 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the 
official policy or position of the Department of Defense or the U.S. Government. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 
Approved for public release; distribution is unlimited. 


13. ABSTRACT (maximum 200 words) 








A plate and frame heat exchanger experimental test stand was developed. Using this test stand a 
performance analysis was conducted. The analysis consisted of evaluating the performance of the heat 
exchanger at varying flow rates and inlet temperatures, to develop an effectiveness-NTU and Log 
Mean Temperature Difference relationships, under steady state operation. The measured heat rates 
were compared to the heat rates provided by the manufacturer and good/bad agreement was found. 
Standard operating procedures for the test stand were developed and implemented. 


15. NUMBER OF 
PAGES * 76 


16. PRICE CODE 


SECURITY CLASSIFI- 18. SECURITY CLASSIFI- ; 19. SECURITY CLASSIFICA- | 20. LIMITATION OF 
CATION OF REPORT CATION OF THIS PAGE TION OF ABSTRACT ABSTRACT 
Unclassified Unclassified Unclassified UL 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 298-102 


SUBJECT TERMS Plate Heat Exchanger. 








— ee 





Approved for public release; distribution is unlimited 


EXPERIMENTAL PERFORMANCE STUDIES OF A PLATE HEAT 
EXCHANGER 


Darren R. Plath 
Lieutenant, United States Navy 
B.A., Macalester College, 1990 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
December, 1996 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOO! 
MONTEREY CA 93943-5101 


ABSTRACT 


A plate and frame heat exchanger experimental test stand was developed. Using this 
test stand a performance analysis was conducted. The analysis consisted of evaluating the 
performance of the heat exchanger at varying flow rates and inlet temperatures, to develop 
an effectiveness-NTU relationship, under steady state operation. The measured heat rates 
were compared to the heat rates provided by the manufacturer and good/bad agreement was 


found. Standard operating procedures for the test stand were developed and implemented. 
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I. INTRODUCTION 


The design of heat transfer equipment is and has been an interesting and important 
subject for some time. In the times of abundant and cheap energy, there was little incentive 
for industry to improve the performance of traditional heat exchangers. With ever increasing 
demands on a limited energy supply, the creation of more efficient and cost effective heat 
exchangers will play a more important role in society. Increased efficiency of heat transfer 
combined with cheaper manufacturing costs are two major factors in the design process of 
heat exchangers. A common design consists of a simple shell and tube configuration. These 
shell and tube heat exchangers are capable of handling a wide range of temperatures and 
pressures, while at the same time requiring a relatively large amount of space. An alternative 
design to the shell and tube heat exchanger that provides a large heat transfer area in a small 
volume is the flat plate and frame heat exchanger. As a general rule, the cost per unit surface 
area decreases with compactness. Flat plate and frame heat exchangers lend themselves to 
high effectiveness designs, thus promoting compactness. As a result, these heat exchangers 
are considered all around the world for many energy conservation and recovery applications. 

Flat plate and frame heat exchangers are constructed of thin plates with relatively low 
thermal resistance. The plates are corrugated and packed close together forming small 
channels, thus taking advantage of the enhanced heat transfer associated with narrow 
passages. These heat exchangers can handle a wide range of fluids, from viscous to non- 
Newtonian. This ability to handle various fluids allows flat plate heat exchangers the 


potential for a wide variety of applications. These heat exchangers are used in the 


automobile industry to control the temperature of paint and to cool automated welding 
machines. There are several other applications including the cooling of bleach solutions and 
boiler blow downs in the pulp and paper industry and process heating and cooling in the food 
industry. In beer breweries, plate and frame heat exchangers are used to cool the beer 
mixture before fermentation. Plate and frame heat exchangers are also used in chemical 
plants, power plants and on maritime shipping vessels. 

The need for a heat exchanger that will resist corrosion and erosion by seawater is of 
prime concern for the Navy. Titanium plate heat exchangers have been found to be highly 
resistant to both conditions even at elevated velocities. In many cases the only required 
maintenance necessitated back flushing the units periodically. With the push to decrease 
manning on Naval vessels, any device that decreases maintenance is advantageous. As a 
result, numerous shell and tube electronic coolers have been replaced with Superchanger 
Navy standard heat exchangers throughout the fleet. 

There are numerous advantages associated with flat plate heat exchangers. They are 
flexible, compact, have low fabrication costs, are easy to clean, afford increased protection 
against erosion and corrosion, provide temperature control with a small temperature 
difference required and have reduced fouling. The disadvantages are potential leakage 
through the gaskets, relatively high pressure drop through the heat exchanger and potentially 
high pumping costs. Another potential down fall of the flat plate heat exchanger is the 
limitation of temperature and pressure. If the temperature or pressure gets too high, there 1s 
a possibility of rupturing the gaskets resulting in leaks. 


The purpose of this project was to produce a system design, component integration, 


and to conduct a performance analysis of a Superchanger plate and frame heat exchanger. 
As part of the performance analysis, the overall heat transfer coefficient and the heat rate 
were derived both experimentally and theoretically and compared with each other. The heat 
transfer coefficient and other performance characteristics were then compared to the 
manufactures value. 

The system consists of a Superchanger UX-056-UJ 23 plate and frame heat 
exchanger manufactured by Tranter inc., two 400 gallon water tanks, a bayonet-style tank 
heater from Cleveland Process Corporation, two Omega FP-6000 paddlewheel flow sensors, 
two single inlet centrifugal pumps powered by Baldor 2-Hp motors and various valves and 
piping. Inlet and outlet temperatures for both the hot and cold sides were measured with 
thermocouples. Numerous runs were conducted with various flow rates for both hot and cold 


sides to evaluate the heat exchanger performance. 





A. 


Il. HEAT EXCHANGER DESCRIPTION 


GENERAL OVERVIEW 

The plate and frame heat exchanger consists of a series of flow passages formed by 
corrugated heat transfer plates, a frame, nozzels, headers and tightening bolts. The 
corrugated plates are held in between the stationary and moveable frames and are compressed 
by the tightening bolts (Figure 2.1). The plates are equipped with elastomeric gaskets and 
have port holes pierced in the corners. When the unit is tightened, the gaskets seal the unit 
and, in conjunction with the nozzels and port holes, allows fluids to flow in alternate 
channels (A and B sides), thus putting the hot fluid and the cold fluid on opposite sides of 
the plate. [Ref. 1]. The dominating function of the headers is to uniformly distribute the 
fluids to the flow passages. 

There are numerous advantages associated with plate heat exchangers. They are easy 
to assemble and disassemble, can be cleaned chemically and mechanically with minimal 
effort and the heat transfer area required can be changed simply with the addition or removal 
of plates. The plate exchanger can be used for multiple services at one time (Figure 2-2), 
floor space required is usually much less than that of a comparable shell and tube exchanger 
and very close temperature approaches can be obtained. In the case of a leak, since the unit 
is open and the gaskets vent to the atmosphere, the fluids will not mix and immediate 
detection is afforded. This can be seen in Figure 2-3. The figure shows three possible leak 
locations in which all the leaks are to the atmosphere. This is especially important when the 


mixing of the fluids could result in damage to machinery. Lube oil coolers are a prime 
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Figure 2.1 Superchanger Flow Diagram 
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Figure 2.2 IsometricFlow Diagram 
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example of this situation where the combining of oil and water is undesirable. Additional 
advantages are low hold up volume, stringent sanitary and hygenic conditions can be met, 
reduced fouling deposits even at low Reynolds numbers and when construction with other 
than carbon steel is required the plate exchanger can be more economical. Probably one of 
the most important advantages 1s the increased heat transfer coefficient that is obtained. This 
increased heat transfer coefficient leads to the requirement for less heat transfer area and in 
turn leads to more compact designs. [Ref. 2] Since the plates are very thin, the conduction 
resistance to heat transfer is minimal, thus providing an additional benefit. 

In spite of all the advantages of the plate and frame heat exchanger there are also 
some disadvantages. The plate exchanger is limited by temperature and pressure. 
Temperatures must be maintained below 175 degrees celsius and pressures above 300 psi are 
unacceptable. These limitations are due to the gasket material seals, plate deformation and 
the frame. If the temperatures become too high, the gasket material and or the glue holding 
the gaskets could melt. This would lead to leakage in the heat exchanger. If the pressure was 
allowed to increase above the stated limits, the gaskets could be blown out of their retaining 
grooves and/or damage to the plates and frame could result. 

Low flow rates, very viscous fluids, and fluids whose viscosity changes radically 
during passage through the unit should be avoided due to the possibility of maldistribution 
of the flow. Another disadvantage is the relatively high pressure drop that occurs in the 
exchanger. As a result, gas to gas service, air coolers and condensing duties can be limited 
as well. [Ref. 2] 


As mentioned earlier, fouling is lower in plate heat exchangers. This is due to the 


highly turbulent flows. This turbulence sustains small particles in suspension. The surface 
of each of the plates is very smooth, despite being corrugated which helps in reducing 
fouling deposits. There is low residence time for crystal formation and outstanding fluid 
dispertion with virtually no low velocity or dead zones. In addition the use of alloys that 
corrode at low rates can be utilized. This practically eliminates any corrosion products from 
fouling the heat exchanger. [Ref. 2] 

The corrugations on the plates form narrow passages for the fluids to pass through, 
and the flow of the fluids is almost always countercurrent. The thin fluid interspace coupled 
with the corrugated plate design induces turbulence that produces extremely high heat 
transfer coefficients. [Ref. 1] The heat transfer coefficient, a, varies as a negative power 
of the flow passage size. A customary expression for the size of a non-circular flow passage 
is the hydraulic diameter, D,, equaling four times the cross-sectional area divided by the 
wetted perimeter. Although a smaller hydraulic diameter increases friction, the benefits of 
compactness on the heat transfer coefficient generally outweigh the detrimental influence of 
small hydraulic diameter on friction. [Ref. 3] 

In addition to the influence of small hydraulic diameter, increases in heat exchanger 
performance can be obtained by any modification of the surface geometry that results in a 
higher heat transfer coefficient at a given flow velocity. One widely accepted modification 
is the use of extended surfaces or fins, so that in addition to providing increased heat transfer 
surface area, the interrupted surface prevents the thickening of the fluid boundary layers. As 
the boundary layer develops, the heat transfer across the boundary layer decreases. 


Therefore, by reducing the boundary layers, the heat transfer is increased. Other methods of 


10 


obtaining increased performance by change of flow surface geometry include the use of 
curved, corrugated, or wavy passages, in which boundary layer separation and turbulence 
(promoters of heat transfer) are induced even at low Reynolds numbers. [Ref. 9] The method 
of corrugations is used in this plate and frame heat exchanger. 

A common descriptor of compact heat exchangers is the area density, which is the 
ratio of heat transfer surface area to heat exchanger volume. A conventional cutoff for 
labeling a heat exchanger as compact is an area density value greater than 700 square meters 
per cubic meter (or 213 square feet per cubic feet). This is not a staunch rule however, as 
many heat exchangers have been grouped into the compact category with lesser area densities 
[Ref. 4]. It will be taken then that compact heat exchangers have a significantly increased 
area to volume ratio compared to that of an orthodox shell and tube unit. The compactness 
of these designs is derived through the need for an extended heat transfer area and in the use 
of plates instead of tubular heat transfer surfaces. 

The corrugations also form contact points between adjacent plates. These contact 
points lend rigidity to the plates, helping them to withstand higher pressures and minimizing 


plate deflections. 


B. PLATE DESCRIPTION 


The plates are manufactured in standard sizes from die-formed sheet metal in 
virtually any material that can be cold worked. Since no welding is used, the weldability of 


the metal is not a concern. They are pressed out of a thin flat sheet of metal. The size, 


1g 


number and arrangement of the plates is contingent upon the duty to be performed. This unit 
has plates with an area of 0.635 square feet per plate. The heat exchanger has 23 plates. 
Two of these plates are in a buffer zone, leaving 21 useable plates. The total heat transfer 
area of the heat exchanger is 13.3 square feet. 

The fluids enter the heat exchanger through inlet nozzles and travel through the port 
holes in the plates. These port holes comprise the header. The individual fluids pass through 
alternating plates in opposite directions. [Ref. 1] Typical plate thicknesses range from .5 to 
Imm. The material used for the plates is chosen on the basis of pressure and corrosion 
conditions. For this plate and frame heat exchanger the plates are made of 316 stainless steel 
and are 0.023622 inches (0.6 mm) thick. 

The corrugations are in a herringbone or chevron pattern (Figure 2-4) and are inclined 
at an angle £ to the direction of flow. The angle § of the unit tested is 60 degrees. Also 
shown in the figure is an alternate washboard plate design. As the angle of inclination 
changes, so does the resulting heat transfer coefficient and the friction factor. The heat 
transfer cofficient increases with B up to an angle of approximately 75 degrees from vertical. 
The friction factor increases by about ten times with an angle change from 30 to 75 degrees. 
[Ref. 5] These corrugations are designed to increase the Nusselt number thus providing 
enhanced fluid heat transfer. The heat transfer enhancement of the corrugations 1s due to 
induced turbulence in the flow. The turbulence eliminates stagnant areas in the liquid flow. 
Since there are no places for the fluids to collect, all of the plate is used, thus increasing the 
effective heat transfer area of the heat exchanger. An additional benefit of the corrugations 


is the added mechanical support to the plate. When an adjacent plate with its chevron pattern 
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Figure 2.4 Plate Details 
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reversed is stacked together, multiple contact points are created. These multiple contact 
points provide support for the thin plates and increase the heat transfer between the hot and 
cold sides of the heat exchanger. [Ref. 1] 

There is a groove around the circumference of the plate which accepts the sealing 
gasket. Due to the nbbed nature of the groove, additional reinforcement is added to the plate. 
The gaskets are a single piece of material, has a molded construction, and is generally bonded 
to the plates with Pliobond 30 adhesive. The gasket material is NBR and is selected based 
on the comapatibility with operating temperatures and the fluids being processed. 

Fluid passage holes are pierced at the corners of each plate. The number and location 
of holes is dependent on the location of the plate in the heat exchanger. Flow directors are 
located at the top and bottom of the heat transfer surface in the port hole areas which evenly 
distribute the fluids. When the heat exchanger is put together, the combination of fluid 
passage holes on the plates comprise a flow path for the fluids. The plates can be arranged 


in a variety of ways depending on the number of passes required. 


Ill. EXPERIMENTAL APPARATUS 
A. PADDLEWHEEL FLOW SENSOR 


The flow sensor chosen for this experimental analysis was the Omega FP-6000 
(Figure 3-1). The fluid flowing by the sensor rotates the paddlewheel. As each blade of the 
paddlewheel passes the magnet in the transducer, it causes a change in the magnetic flux field 
running through the coil and core. When this flux field changes, an alternating current (AC) 
voltage, approximating a sine wave, 1s induced. The frequency of the output voltage of the 
coil is directly proportional to the linear velocity of the fluid in the pipe; therefore, it 
represents the amount of fluid passing the flow sensor in a particular increment of time. A 
complete cycle occurs each time one of the paddlewheel blades passes the coil; therefore, 
four complete cycles are generated for each paddlewheel rotation. The FP-6000 senses local 
velocity at an insertion depth of 10% of inside diameter, therefore, the velocity that it detects 
iS approximately 95% of the mean velocity in the pipe. The flow sensor is constructed of 
brass and is completely self contained with out any user serviceable parts. [Ref. 6] 

Installation of the flow sensor started with the calculation of the “H” dimension 
(Figure 3-2). The following equation was used: 

H = 5.95 - wall thickness - 10% of inner diameter (3.1) 
The ‘“H” dimension was found to be 5.55 inches. Next the sensor was located in an 
undisturbed section of pipe at the twelve o’clock position just prior to the entrance to the 
plate heat exchanger. The lower hex nuts were adjusted to allow the moving of the flow 


sensor to the calculated ““H” dimension. The jam nuts were tightened securing the sensor. 
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Figure 3.1 Cross Section of the FP-6000 
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Figure 3.2 FP-6000 Dimensions for Calculating “H” Dimension 
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Figure 3.3 Downstream Alignment using 10" Metal Rod 
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The sensor was then inserted in the saddle fitting with the alignment hole on the flange 
pointing in the down stream direction of the flow. A ten inch rod was provided with the flow 
sensor to make minor adjustments to the alignment of the sensor (Figure 3-3). The wire 
leads from the flow sensors were fitted with BNC plugs. A digital frequency reader was used 


to measure the flow rates. 


B. |BAYONET-STYLE TANK HEATER 


In order to achieve a sufficient temperature difference, an imersion style tank heater 
was used. A Cleveland Process Corporation Bayonet-style imersion tank heater was 
selected. The heater is constructed of 304 stainless steel, has three elements, is three phase 
480 volts , and 27,000 watts. The heater has an auto shut off device incase a low water 
situation arises. A bracket was manufactured to hang the tank heater in the middle of the 
tank. The bayonet-style tank heater is 51 inches in length. The hot zone extends 42 inches 
from the bottom. 


The required heat rate of the heater was determined using the following equation: 
Q=m CAGAT (3.2) 
The mass of the water was calculated by multiplying the volume of the tank with the density 
of water at room temperature. Room temperature was taken to be 22 degrees Celsius. The 


temperature difference used in the previous equation is, the difference in temperature 


between the water leaving the tank and that of the water returning to the tank. A maximum 
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temperature difference of 10 to 15 degrees is expected. With AT set at the upper limit of 15 
degrees, the resulting heat rate Q, is 113.452 MJ. While this requirement is larger than the 
heater selected, it is assumed that over a relatively short duration, the temperature will remain 
fairly constant in the tank. Since the tank is rather large and the heater is located between 
the returning water and the out going water, this assumption should be valid. It is noted that 
the majority of runs will be conducted with the delta T being less than 15 degrees. In 


addition, as the flow rate increases, the time duration of constant temperature will decrease. 


C; SELF-PRIMING CENTRIFUGAL PUMP 


Two single inlet centrifugal pumps powered by Baldor 2-Hp motors were chosen for 
this experimental performance analysis. The pumps are self priming and are designed for 
continuous transfer pumping of liquid from tanks or sumps. Each pump has twelve feet of 
suction lift. The motor operates at 3450 RPM, 115 volts and 60 Hz. 

A preliminary system headloss calculation was conducted using the Darcy-Wiesbach 


equation: 


ae (3.3) 


Minor losses were also accounted for. K values for 90 degree elbows were taken to be 1.4, 
for valves k equaled 0.26. [Ref. 7] The pressure drop through the heat exchanger used was 


6 psi (from manufacturer). Measurements of pipe lengths and number of elbows and valves 
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were taken from the preliminary system diagram. The piping used was schedule 40 PVC. 
The pipe is considered smooth and a friction factor of 0.015 was used. The head loss 
calculated was 32.64 ft of water. 

In planning the performance analysis, it was decided that a wide range of flow rates 
would be desirable. Cross referencing several vendor catalogs provided us with several 
options to choose from. The pump chosen, is capable of providing 115 gpm at 35 ft of head. 


This is more than sufficient for this experimental performace analysis. 


D. SYSTEM CONSTRUCTION 

The experimental apparatus started with a preliminary piping diagram. From this 
diagram measurements were taken to calculate the length of pipe and the number of valves 
required. Piping, valves and assorted other items were then purchased. The two water tanks 
were cleaned and positioned in the laboratory. Once all items were received construction 
began. 

The first step involved drilling holes in the floor to secure one of the pumps in place. 
The inlet piping was then connected. Each threaded connection was wrapped with teflon 
tape, to protect against the possibility of leaks. Non-threaded connections were sanded, 
primed and then glued together. The next step was to hook up the discharge of the pump to 
the inlet of the heat exchanger. The paddlewheel flow sensor and thermocouple bead were 
inserted in this section. A recirculation line was also installed. Next a line was connected 
from the discharge of the heat exchanger to the tank. This completed one side of the heat 


exchanger setup. The other side was identical. The paddlewheel flow sensors were then 
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positioned according to the operating instructions. Before the thermal couple beads were 
inserted to the system, they were calibrated. The thermocouple calibration is discussed in 


Appendix B. 
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IV. SYSTEM ALIGNMENT AND OPERATION 


A. INTRODUCTION 
The objective of this chapter is to discuss the basic philosophy and decision making 

process used in developing the Standard Operational Procedures for the plate and frame heat 
exchanger system. The primary goal in the development of the standard operational 
procedures was to produce a detailed step by step procedure which can be followed by an 
operator with a general knowledge of the overall system. The major criteria are both safety 
of the operator and to prevent damage to equipment. Two procedures have been developed 
to help accomplish the above goals. They are listed below: 

* Heat Exchanger System Recirculation Procedure (HESRP) 

* Heat Exchanger System Operational Procedure (HESOP) 


The above written procedures are contained in Appendix A. 


B. SYSTEM ALIGNMENT 

As mentioned earlier the experimental apparatus consists of two water tanks, two 
centrifugal pumps, two paddlewheel flow sensors and a plate and frame heat exchanger. A 
schematic diagram of the system is shown in figure 4-1. The hot and cold sides of the system 
are exactly the same with the exception of valve numbering. HW and CW will be used to 
denote the hot and cold side valves respectively. For this reason, only one side will be 
dicussed. 


The flow path of the water will be described starting at the water tank. Water flows 
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Figure 4.1 System Schematic Diagram 
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through the pump suction valve HW-1 into the centrifugal pump. From the pump it flows 
through the discharge valve HW-2. After this valve the line branches in two. One branch 
is a recirculating line back to the tank. This recirculation line will allow for increased 
mixing of the fluid in the tank prior to each test run. This mixing will help allow the 
assumption of a constant temperature throughout the tank to be valid. The other line 
continues toward the heat exchanger. The water passes through the paddlewheel flow sensor 
and into the heat exchanger unit. After traversing the heat exchanger the fluid passes through 
valve HW-4 and then back to the tank. Valves HW-4 and CW-4 will serve as throttle valves, 
permitting the capability to increase or decrease the flow rates as necessary. Temperature 
measurements are taken prior to entering and upon exiting of the heat exchanger as well as 


in both tanks. 


Cc: OPERATING PROCEDURES 

The procedure starts with a safety space walk through, system verification and system 
alignment. During the visual inspection of the heat exchanger lab space, emphasis is placed 
on removing water and any debris from the area, as well as, checking for leaks and damaged 
equipment. After completion of the safety walk through the operator commences with 
placing the system in recirculation mode. The operation of the system is quite simple, 
however the correct valves must be open for the water to flow through the heat exchanger. 
Damage to the heat exchanger could result from improper operation. 

The first step is to turn on the heater. The water must be heated to the desired starting 


temperature. Before allowing any fluid to flow the operator must ensure the plate pack 
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length is correct. The maximum length is 2.9 inches and the minimum length 1s 2.7 inches. 


Figure 4-2 shows the plate pack thickness as the dimension “A”. Once the water has reached 





Figure 4.2 Plate pack thickness 


the required temperature, the pumps can be aligned for operation. The operator will first 
place the system into a recirculation mode in accordance with procedure (HESP) in 
Appendix A. The following paragraph gives a brief explanation of this procedure. 
Ensure that both HW-1 and CW-1 are open. Second, close valves HW-2 and CW-2 
and open valves HW-3 and CW-3. Next check to see that valves HW-4 and CW-4 are 
closed. As a reminder both drain valves HW-5 and CW-5 are to remain closed at all times 
during the experimental analysis. The next step is to start both pumps. Once both pumps 
have been started, slowly open the discharge valves HW-2 and CW-2. This will allow the 


fluid to flow through out the sytem in a recirculating mode. 
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When ready to procede with the experimental analysis, the operator will place the 
system in operation in accordance with procedure (HESOP) Appendix A. Slowly open 
valves HW-4 and CW-4. Open the cold side first. This will ensure that cooling water is 
flowing through the heat exchanger prior to the hot fluid. For this experimental analysis, the 
cold fluid is room temperature water. Since the heat exchanger will be at room temperature 
this procedure will reduce the possibility of any thermal shock or over heating of the heat 
exchanger. Now that the system is in operation, use valves HW-4 and CW-4 to adjust the 
flow rate through the heat exchanger. 

It is important to monitor the system for any possible leaks that may occur. In the 
case of a leak, the operator needs to identify which side of the system it has occured on. If 
a leak occurs, slowly shut both pump discharge valves HW-2 and CW-2, then stop the 


pumps. Ensure corrective measures have been completed before operating the system again. 


Mi 





V. DATA COLLECTION AND ANALYSIS 


re APPROACH 

In conducting the performance analysis of the plate and frame heat exchanger, the 
Effectiveness - NTU and Log Mean Temperature difference methods were utilized. The 
major requirement for these methods is the determination of a heat exchanger overall heat 
transfer coefficient. This includes the computing of other parameters such as the hot and 
cold side convection coefficients, wall conduction resistance, and fouling resistance. The 
overall heat transfer coefficient was found both from theoretical and experimental data. 
Since low fouling factors are an inherent advantage of plate and frame heat exchangers, the 
calculation of a fouling resistance will be neglected, however the negative impact of fouling 


will be discussed. 


B. OVERALL HEAT TRANSFER COEFFICIENT 

1. Newton’s Law of Cooling 

The total thermal resistance to heat transfer between two fluids separated by a heat 
exchanger structure is termed the overall heat transfer coefficient, U. The determination of 
U is essential to any heat exchanger analysis. In an expression analogous to Newton’s Law 
of Cooling, the heat transfer between two fluids separated by one or more thermal 


resistances, R,, 1S: 


Q = =— = UA AT G1 


Py) 


The result apparent from the above expression, that 1/).R, = UA,, is applicable to 


clean, unfinned surfaces only, but can be used as the basis for determining a heat exchanger’s 
overall heat transfer coefficient with additional factors considered such as the effects of wall 
conduction and the presence of fouling on both hot and cold water sides. As mentioned 
earlier, the effects of fouling will be neglected in the calculations. [Ref. 8] 

746 Convection Coefficients 

The water side convection coefficients, a, and @ , are computed in the following 
manner as recommended in Hewwit’s Process Heat Transfer. [Ref. 5] 

(1) Mass Flow Rate. The mass flow rate per passage is found from the 


following equation: 


m = (5.2) 





where, M is the total mass flow rate of the system. 
(2) Cross Sectional Area. The cross sectional area of each passage is given 


by: 


S = bW (257) 


where 6 is the gap between plates and W is the width of the plate. 
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(3) Water Velocity. The water velocity across the plate is given by: 


“= — (5.4) 


(4) Reynolds Number. The Reynolds number, Re, can be computed using the 
following relation: 


uD, 
Re = ——— (5.5) 
yl 


(5) Nusselt Number. The approximation for the Nusselt Number, Nu, in the 


turbulent region is 


Nu = 0.4 (Pr)** (Re)? (5.6) 


(6) Convection Coefficient. The convection coefficient can now be found 


using the following relation: 





a= (5.7) 


where &, is the thermal conductivity of the fluid. 

(7) Sequence of Computations. Computations of the convection coefficients 
for each side of the system’s heat exchanger were completed, and are shown as part of 
Appendix D. 

5: Fouling Resistance 
During normal heat exchanger operation, surfaces are often subject to fouling by fluid 


impurities, rust or scale formation, and other reactions between the fluid and the wall 
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material. Depositions such as these can greatly increase the resistance to heat transfer 
between the fluids. This effect can be treated by introducing an additional thermal resistance, 
termed the fouling factor, Ry the value of which depends on the operating temperature, fluid 
velocity, and length of service of the heat exchanger. The heat exchanger used is brand new 
and the amount of fouling expected is extremely small. Since the values for the fouling 
factor are very small compared to the convection coefficients and the thermal resistance of 
the plates, they can be neglected. 

4. Plate Conduction 

Derivations of the plate conduction resistance, R,,,, by Incropera and DeWitt [Ref.7] 


produce the following expression for heat exchanger plates (modeled as plane wall): 


(5.8) 


where f¢ is the plate thickness and & is the thermal conductivity of the plate. 
5. Summation of Elements 
The summation of thermal resistances including effects of hot and cold side convec- 


tion, and plate conduction yields the following: 


1-1, 2 (5.9) 
UO, OL, 


The area is not included in equation 5.9 since the area is the same for both sides of the plate. 
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c. THE EFFECTIVENESS - NTU METHOD 

1. Applicability 

The Effectiveness - NTU method utilizes the overall heat transfer coefficient to 
determine heat exchanger heat transfer performance characteristics. In the Effectiveness - 
NTU approach, with knowledge of the heat exchanger type, size, and fluid flow rate, the 


NTU, maximum heat transfer rate, 9 __, and effectiveness, €, the actual heat transfer rate 
max 


can be determined. 

Ze Maximum Possible Heat Transfer Rate 

The maximum possible heat transfer rate for a given heat exchanger is achieved 
when, by definition, the entering cold fluid temperature reaches the maximum possible 
temperature difference, 1.e., that of the entering hot fluid temperature. Using capital T for 
the hot stream and little t for the cold stream, then T1 and T2 correspond to the entrance and 
exit temperatures of the hot stream respectively. Alternately, t3 and t4 correspond to the 
entrance and exit temperatures of the cold 
stream respectively. Thus the entering cold fluid would rise in temperature by the quantity 


(t, - t,). The heat transfer rate can be found using the following relationship: 
OBnC(t = t,) (5.10) 
where the term mc, is defined as the heat capacity rate, C, and will be based on the mass 


flow rate, m , and specific heat, c >? of the fluid. The maximum heat transfer rate, O ax ; 


£8) 


is then: 
O ax ~ C min (Fy a ts) (5.11) 


whereC,,, is equal to eitherC,,, or C,,,,, whichever is smaller. [Ref. 8] 


cold? 


3. Effectiveness 
The effectiveness, €, is defined as the ratio of the actual heat transfer rate for a heat 


exchanger to the maximum possible heat transfer rate. 


€ = (5.12) 


aoe 
One 

The effectiveness, which is dimensionless, must be in the range 0 < € < 1 .[Ref. 5] 
It is useful because, if e, T,, and ¢, are known, the actual heat transfer rate may be 


determined from the expression: 


Q=€0O. =€C..(T, - t;) (5.13) 


4. Number of Transfer Units 
The number of transfer units, NTU, as explained by Incropera and Dewitt [Ref. 8], 
is a dimensionless parameter that is widely used for heat exchanger analysis and is de- 


fined as:[Ref.5] 





NT a (5.14) 
C 
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=p Effectiveness - NTU Relations 
The Effectiveness - NTU relationship for the type of heat exchanger being ana- 
lyzed, i.e. a single pass counter flow plate heat exchanger with both fluids unmixed is: 


expan. © )NTU | - 1 
E = (5:15) 
expul(ieese ) NIU | 


where: 


_—— (5.16) 





D. THE LOG MEAN TEMPERATURE DIFFERENCE METHOD 


The log mean temperature difference method makes use of parameters such as the 


inlet and outlet temperatures and the mass flow rate of the fluids. With inlet and outlet 


temperatures and mass flow rates known, the heat rate can be obtained by applying overall 


energy balances to the hot and cold fluids. If Q is the total rate of heat transfer between the 


hot and cold fluids and there is negligible heat transfer between the exchanger and its 


surroundings, as well as negligible potential and kinetic energy changes, application of an 


energy balance gives: 


Onan aCr. (Dal, (5.17) 


and 


ORs Cpl... -T...)) (5.176) 


where the temperatures appearing in the expressions refer to mean fluid temperatures at the 
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designated locations. [Ref. 8] 
Another way of looking at the relationship between the total heat transfer rate Q and 


the temperature difference, AT, is: 


Q = UAFoT, (5.18) 


where AT,,, is an appropriate mean temperature difference and F is a correction factor to 
account for the non-counter flow configuration. The correction factor, F, is a function of the 
ratio of mass flow rates, temperature differences and heat exchanger configuration. The 
value for F was taken to be 0.96 for this particular heat exchanger. [Ref. 5] For plate heat 
exchangers with countercurrent flow and both fluids unmixed, the form of AT, is as 


follows: 


= 20) , Peat ) 


AT (5.19) 
1 alee 
n—_--————- 
i le 


From equation 5.18 the overall heat transfer coefficient can be found using equations 5.17 


and 5.19. 
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VI. DISCUSSION OF RESULTS 


A. INTRODUCTION 

From the data collected three important heat exchanger parameters were calculated 
and compared in this experimental study. They are the overall heat transfer coefficient, heat 
exchanger effectiveness and the heat rate. All values were first calculated from correlations 
recommended by Hewitt in Process Heat Transfer [Ref. 5] and then from the measured 
experimental data. A detailed discussion of these parameters follows. 

B. OVERALL HEAT TRANSFER COEFFICIENT 

With knowledge of the mass flow rate and the inlet temperatures of the hot and cold 
fluids, the Reynolds number, Nusselt number and finally the convection coefficient for both 
the hot and cold sides of the heat exchanger were calculated from available correlations. The 
convection coefficients along with the plate conduction comprise the overall heat transfer 
coefficient. This gives us the expected value of the overall heat transfer coefficient. 

The overall heat transfer coefficient was also found by experimentation using the log 
mean temperature difference method. Using the simple energy balance of equations 5.17a 
and 5.17b the heat rate for the hot and cold sides were computed. With the known inlet and 
outlet temperatures along with the heat rate obtained, equation 5.18 can be used to obtain 
UA. The value of F was given earlier to be 0.96. The log mean temperature difference was 
computed using equation 5.19. 

The experimental values for the overall heat transfer coefficient were approximately 


15 -20% lower than the expected values. This can be attributed to using an a correlation for 


a7 


the Nusselt equation based on a standard plate exchanger and from the extraction of the 
experimental heat transfer coefficient from the energy balance. The Nusselt correlation could 
produce a slightly higher value while the energy balance yields a lower value. The energy 
balance will be discussed in the heat rate section. 

All calculations were performed in a spreadsheet and are included as Appendix D. 

C. HEAT RATE 

The heat rate was calculated in two different ways. The first utilized the Log Mean 
Temperature Difference method, and the second finds the heat rate using the simple energy 
balance. Using the value of the overall heat transfer coefficient calculated earlier, and 
equation 5.18, with the area provided by the manufacture and the correction factor discussed 
earlier, the expected value of the heat transfer rate can be obtained. 

The second method used to acquire the heat rate was the energy balance equations 
of 5.17a and 5.17b. In principle, the heat transfer from one fluid must equal the heat transfer 
to the other fluid. The heat rates for the cold and hot fluids were compared and found to be 
well within the estimated uncertainty limits. The uncertainty ranged from a low of 4.62% 
to a high of 28.38%. The high value of uncertainty is attributed to cases of large disparity 
between the flow rates of the hot and cold fluid streams. For the most part the values were 
within 10% of each other. 

The experimental heat rate values were less than the expected estimates. A small 
portion of the difference can be accounted for by the assumption that all the heat is 
transferred between the fluids. This assumption 1s not entirely true since some heat is lost 


through the heat exchanger to the surroundings. A standard fin analysis was conducted on 


38 


the heat exchanger and 1t was discovered that the losses were on the order of 100 watts. A 
second assumption used was that the fluids are evenly distributed between the channels. It 
was also assumed that the heat exchanger was in steady state operation. 
D. HEAT EXCHANGER EFFECTIVENESS 

As seen in figure 6-1, the effectiveness-NTU relationship in equation 5.15 is used to 
produce the solid and dotted line curves for increasing Cr with the lower bound having a Cr 
value of .999 and the upper bound value of 0. The experimental data fits within these bounds 
very nicely. The dots correspond to values of Cr between 0.97 and 0.99 while the triangles 
correspond to various Cr values ranging from 0.17 to 0.80. It can be seen from figure 6-1 
that the theoretical values obtained from the effectiveness NTU relation correlate well to the 
experimental data. Figure 6.2 shows the same Effectiveness-NTU plot with the experimental 
uncertainty superimposed as error bars. The larger uncertainties are due to large disparities 


in the flow rates. 
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Figure 6.1 NTU vs. Effectiveness 
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Vil. CONCLUSIONS AND RECOMMENDATIONS 


The system design, component integration and performance analysis of the 
Superchanger plate and frame heat exchanger has been completed. Various performance 
characteristics were calculated from recommended correlations and measured experimental 
data. The overall heat transfer coefficient, the heat rate and the effectiveness of the heat 
exchanger were deduced from experimental data and compared to values calculated from 
correlations. The comparisons have been described earlier. From the discussion of the 
results it has been determined that the correlations can be used to predict the results of the 
heat exchanger satisfactorily. The expected values will be slightly higher than the 
experimental measurements. When the manufacture’s data was compared to the 
experimental data, it was found that the manufacturer’s heat rate and heat transfer coefficient 
were also higher than the experimental results. 

Additional data runs should be conducted with ethylene glycol as one of the operating 
fluids. This would allow for greater inlet AT’s and continuous cooling of the cold tank to 
maintain steadier temperatures. With the ability to continuously cool the cold tank, the 
length of each run would then be increased. With longer runs the insurance of steady state 


can be met. 
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APPENDIX A. STANDARD OPERATIONAL PROCEDURES FOR 
ALIGNMENT AND OPERATION OF THE SUPERCHANGER 
UX-056 PLATE AND FRAME HEAT EXCHANGER 


HEAT EXCHANGER SYSTEM RECIRCULATION PROCUDURE 
(HESRP) 


PROCEDURE 


SYSTEM ALIGNMENT FOR RECIRCULATION 


1. Conduct visual inspection of heat exchanger system and verify the following: 
a. Ensure the deck are free of water and dry. 
b. Inspect all piping runs and asccessories for loose connections, damage, or leaks. 
c. Ensure all valve handles are installed and valve labels are in place. 
d. Check water level in both tanks. Ensure enough water is present to support heat 
exchanger operations and ensure water level in hot water tank is above water level 


cutoff sensor on tank heater. 


2. Ensure the following valves are in the fully open position: 


a. Hot side pump suction valve HW-1 
b. Cold side pump suction valve CW-] 
c. Hot side recirculation valve HW-3 
d. Cold side recirculation valve CW-3 
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3. Ensure the following valves are in the fully closed position: 


a. Hot side pump discharge valve HW-2 
b. Cold side pump discharge valve CW-2 
c. Hot side heat exchanger discharge valve HW-4 
d. Cold side heat exchanger discharge valve CW-4 
e. Hot water tank drain valve HW-5 
f. Cold water tank drain valve CW-5 


NOTE: Both drain valves must remain shut during operation of pumps and heat exchanger. 


4. Start the cold side pump. Once pump has been started, slowly open pump discharge 
valve CW-2 

5. Start the hot side pump. Once pump has been started, slowly open pump discharge 
valve HW-2 


6. Monitor system for possible leaks. 
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HEAT EXCHANGER SYSTEM OPERATIONAL PROCUDURE 
(HESOP) 


PROCEDURE 


CAUTION 
Check and maintain required water level in hot side water tank. Low water level could 
result in over heating of tank heater. Prior to operation of system monitor temperature of 


hot water tank to ensure no damage occurs to tank or tank heater. 


SYSTEM ALIGNMENT FOR OPERATION 

1. Align and place system in recirculation (HESRP) 

2. Crack open cold side heat exchanger discharge valve CW-4. 

3. Slowly close cold side recirculation valve CW-3. 

4. Crack open hot side heat exchanger discharge valve HW-4. 

5. Slowly close hot side recirculation valve HW-3. 

6. Monitor system for possible leaks. 

7. If leak occurs on hot side, stop hot side pump and close pump discharge valve HW-2. 
8. If leak occurs on cold side, complete step 7. Then stop cold side pump and close pump 


discharge valve CW-2. 
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APPENDIX B. THERMOCOUPLE CALIBRATION 


In the field of heat transfer it is vital that temperature measurements be accurate. To 
ensure accurate measurements are made, calibration of temperature reading devices are needed. 
This appendix outlines the procedure used in the calibration of the T-type thermocouples used 
in the experimental analysis of the Superchanger UX-056-UJ plate and frame heat exchanger. 

The thermocouple wire was cut to the specified length. Beads were then made using a 
Dynatech Corporation Thermocouple Welder. The thermocouple beads were placed in a Keithly 
740 System Scanning Thermometer. All beads were then immersed in a Rosemount Engineering 
Model 913A ethyl glycol bath. An associated Rosemont Engineering Model 920A commutating 
bridge provided the reference temperature as measured by a precision temperature probe. A four 
wire bridge electrical circuit was used to keep the lead resistance as low as possible. 

The temperature of the ethyl glycol bath was slowly increased in small increments from 
roughly 15°C to 75° C. This covers the experimental temperature range of interest. When the 
temperature of the bath steadied out at each incremental change, the corresponding resistance 
from the commutating bridge and the reading from each channel on the thermometer scanner was 
recorded. Table B-1 lists all of the data recorded. The channel numbers correspond to individual 
thermocouple beads. The rough temp was used to make the incremental temperature changes. 

The calibration of the T-type thermocouples yielded results within 0.1° C of the reference 
temperature up to a temperature of 50°C. Above that the readings were within 0.3 C of the 


reference temperature. 
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APPENDIX C: UNCERTAINTY ANALYSIS 


A. Velocity. 


The equation for velocity is 


m 
“a4 -e—_— 








(C.1) 
p 
where 
Ciara (ED 
Noto 
and 
M = PQ (C3) 
The uncertainty becomes 
Ou ) 
U = — U C.4 
3 (Ss 1 3 
Of 10) 
oye (C.5) 
u O 


B. Reynold’s Number 


The Reynold’s number equation 1s 
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Ress ee (C.6) 
u 
The uncertainty becomes 
7 ORe : 
W = ( u,] (ED 
U U 
a - = (E28) 
C. Nusselt Number 
The Nusselt number equation is 
Nu = 0.4 (Pr)°* (Re)?™ (C.9) 
The uncertainty becomes 
7 ONu Z 
Uru = || ap, Ure (C.10) 
Uy. 0.64 U 
ees (C.11) 
Nu Re 
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D. Heat Rate 


The equation for the heat rate is 


Operon T, ) (C.12) 


The uncertainty becomes 





(En3) 
(C.14) 
E. Effectiveness 
The equation for the effectiveness is 
e = £ (C.15) 
Cae 
(C.16) 





and 


38) 
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APPENDIX D. EXPERIMENTAL DATA 


This appendix contains all pertinent data that was collected and the spread sheets that 


contain all the calculated parameters. 


Se 








43.8) 35.5) 





| 


4 























6576-04 


























14 









































0. 085479 





11.9792 




































































































































































1 20a) 14. o740 0,000946 , 0.940264 | di 5000756 "995, 
2 | 43.0| 34.3; 24.0/31.8) 23 | 22 | 17.2201 | 0.001086 6.70E-04 | 1.081323} 0.098302 | 16.4714 | 0.001039 | 995.503 
3 | 41.6| 35.8| 26.7/33.0! 31 | 30 | 23.2097 | 0.001464 6.70E-04 | 1.457434 | 0.132494 | 22.4610 | 0.001417 | 995.469 
4 | 42.5| 36.9] 29.1/34.2) 33 | 35 | 24.7071 | 0.001559) 6.57E-04 | 1.551434| 0.141039 | 26.2045 }- 0.001653 | 995.437 
5 | 43.8 | 39.2} 31.6/36.4| 41 | 41 | 30.6967 | 0.001937| 995.26 | 6.35E-04 | 1.927478] 0.175225 | 30.6967 | 0.001937 | 995.395 
6 | 43.2} 40.0] 34.2/37.8) 46 | 45 | 34.4402 | 0.002173] 995.26 | 6.33E-04 | 2.162532] 0.196594 | 33.6915 | 0.002126 | 995.359 
7 =| 44.2) 40.3| 35.0'39.0] 26 | 26 | 19.4662 | 0.001228| 995.25 | 6.26E-04 | 1.222286 0.111117 | 19.4662 | 0.001228 | 995.341 
8 | 42.8] 35.3/ 25.6/33.5| 22 | 22 | 16.4714 | 0.001039| 995.3 | 6.65E-04 | 1.034302] 0.094027 | 16.4714 | 0.001039 | 995.474 
9 | 43.3/ 38.0| 27.6 | 33.2; 53 | 54 | 39.6811 | 0.002503) 995.28 | 6.45E-04 | 2.491655| 0.226514 | 40.4298 | 0.002551 | 995.459 
10 | 42.6/ 39.4| 31.4/35.1) 88 | 84 | 65.8856 | 0.004157 995.27 | 6.41E-04 | 4.137062| 0.376097 | 62.8908 | 0.003968 | 995.408 
11 | 44.9] 38.1] 27.9; 35 | 30 | 30 | 22.4610 | 0.001417) 995.26 | 6.35E-04 | 1.410349| 0.128214 | 22.4610 | 0.001417 | 995.438 
12 | 46.3] 40.1| 29.5) 36 | 40 | 42 | 29.9480 | 0.001889| 995.23 | 6.15E-04 | 1.880409/ 0.170946 | 31.4454 | coeese 995.417 
13 | 44.8) 41.6| 35.2} 39 | 59 | 59 | 44.1733 | 0.002787| 995.23 | 6.15E-04 | 2.773603) 0.252146 | 44.1733 | 0.002787 | 995.343 
14 | 42.0) 32.9| 23.9) 33 | 14 | 15 | 10.4818 |0.000661| 995.33 | 6.86E-04 | 0.658211| 0.059837 | 11.2305 | 0.000709 | 995.498 
15 | 42.0| 37.0| 245| 35 | 31 | 15 | 23.2097|0.001464| 995.3 | 6.59E-04 | 1.457413| 0.132492 | 11.2305 | 0.000709 | 995.468 
16 | 42.9] 39.4| 25.1| 37 | 48 | 15 | 35.9376 | 0.002267) 995.27 | 6.39E-04 | 2.256573| 0.205143 | 11.2305 | 0.000709 | 995.449 
17 | 43.0| 40.0| 24.1) 36 | 70 | 17 | 52.4090 | 0.003306] 995.26 | 6.35E-04 | 3.290815| 0.299165 | 12.7279. 0.000803 | 995.462 
18 | 41.7| 39.8| 25.6! 36 | 93 | 17 | 69.6291 | 0.004393| 995.27 | 6.44E-04 | 4.372142! 0.397467 | 12.7279 | 0.000803 | 995.451 
19 | 41.6} 32.4| 27.4/31.1| 13 | 31 | 09.7331 | 0.000614] 995.34 | 6.92E-04 | 0.611201] 0.055564 | 23.2097 | onto 995.479 
20 | 41.9| 32.7| 28.7|31.1| 13 | 49 | 09.7331 | 0.000614] 995.34 | 6.88E-04 | 0.611197) 0.055563 | 36.6863 | 0.002315 | 995.468 
21 | 42.0| 32.9] 29.7/31.2| 13 | 77 | 09.7331 |0.000614| 995.33 | 6.86E-04 | 0.611196) 0.055563 | 57.6499 | 0.003637 | 995.458 
22 | 45.6 | 36.9| 23.8 | 32.9] 31 | 32 | 23.2097 |0.001464| 995.27 | 6.38E-04 | 1.457368/ 0.132488 | 23.9584 | 0.001512 | 995.495 
23 | 45.4| 39.3/ 24.6|34.7| 49 | 32 | 36.6863 | 0.002315] 995.25 | 6.25E-04 | 2.303536] 0.209412 | 23.9584 | 0.001512 | 995.472 
24 | 44.3| 40.4| 27.1/36.1] 68 | 32 | 50.9116 | 0.003212| 995.25 | 6.25E-04 | 3.196743 0.290613 | 23.9584 | 0.001512 | 995.438 
25 | 42.8 40.3| 30.2|37.1| 88 | 30 | 65.8856 | 0.004157) 995.26 | 6.34E-04 | 4.137021| 0.376093 | 22.4610 | 0.001417 | 995.401 
26 | 30.1| 23.5| 21.4] 22 | 10 | 85 | 07.4870 | 0.000472| 995.52 | 8.60E-04 | 0.470240} 0.042749 | 63.6395 | 0.004015 | 995.612 
27 | 35.5| 32 | 245} 28 | 66 | 63 | 49.4142 | 0.003118] 995.4 | 7.40E-04 | 3.103199| 0.282109 | 47.1681 | 0.002976 | 995.532 
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0.47956 





0.18695 | 0.17882 
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0.37390 | 0.36577 
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0.71528 0.68277 
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0.56897 | 0.13818 


0.75592 | 0.13818 




















8.35E-04 0.06840 
8.39E-04| 1.03451 | 0.09405 
8.04E-04| 1.41064 | 0.12824 
7.73E-04| 1.64570 | 0.14961, 
7.36E-04 | 1.92774 | 0.17525 
7.06E-04 | 2.11573 | 0.19234 
6.92E-04 | 1.22240 | 0.11113 
8.09E-04| 1.03448 | 0.09404 
7.94E-04| 2.53913 | 0.23083 | 
7.47E-04| 3.94956 | 0.35905 | 
7.74E-04 | 1.41060 | 0.12824 
7.55E-04) 1.97480 | 0.17963 
6.94E-04 | 2.77392 | 0.25217 
8.33E-04 0.70534 | 0.06412 
8.03E-04 | 0.70532 | 0.06412 
7.85E-04 | 0.70531 | 0.06412 
7.98E-04 | 0.79936 | 0.07267 
7.86E-04 | 0.79935 | 0.07267 
8.14E-04 | 1.45768 | 0.13252 
8.03E-04 | 2.30405 | 0.20946 
7.93E-04 | 3.62061 | 0.32915 
-8.31E-04| 1.50473 | 0.13679 | 
8.07E-04| 1.50469 | 0.13679 | 
7.74E-04 | 1.50464 | 0.13679 
7.41E-04| 1.41055 | 0.12823 
9.66E-04 | 3.99740 | 0.36340 
8.70E-04 | 2.96254 | 0.26932 


0.10567 | 0.39828 
0.10567 | 0.62587 





0.25197 | 0.26010 
0.39828 | 0.26010 


0.55272 0.26010 





0.10567 | 0.25197 





0.71528 | 0.24385 
0.08128 | 0.69090 





0.53646 | 0.51208 


| 1,279.868 
| 1,443.568 | 1,103.514 
1,947.577 | 1,570.727 
2.113.818 | 1,904.789 
2.717.950 | 2,343.932. 
| 3,055.160 | 2,679.752_ 
1,748.020 | 1,579.631 | 
1,391.600 | 1,144.368 | 
—3,457.289 | 2,861.184 
5,778.732 | 4,728.059 
1,988.744 | 1,630.942 
2,737.163 | 2,339.448 
| 4,037.316 | 3,677.413_ 
0,858.205 | 0,757.427 
1,978.059 | 0,786.217 
3,161.011 | 0,804.231 | 
| 4,640.402 | 0,896.859 
-6,078.126 | 0,909.509 
-0,789.816 | 1,601.984 
| 0,794.539 | 2,568.309 | 
-0,796.905 | 4,084.240 
—-2,045.354 | 1,621.224 
3,300.504 | 1,668.168 


4,580.292 | 1.739.671. 
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76.4617 | 70.9367 
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97.0709 | 97.2517 
112.1211 | 108.4569 
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6,940.330 3.95E “05 | 
8.485.122 | 3.95E-05| 3,749.660 
10,515.069 | 3.95E-05| 4,447.027 
~11,632.796 | 3.95E-05| 4,726.170 
12,973.118 | 3.95E-05 | 5,231.607 
13,900.676 | 3.95E-05| 5,534.315 
9,826.594 | 3.95E-05| 4,155.556 
8,610.857 | 3.95E-05 | 3,731.749 
15,400.988 | 3.95E-05 5,956.875 
20,540.749 | 3.95E-05| 7,456.667 
10,565.045 | 3.95E-05| 4,440.174 
13,176.120 | 3.95E-05| 5,239.530 
16,595.411 | 3.95E-05| 6,305.111 
6,664.258 | 3.95E-05/ 2,942.921 
6,752.348 | 3.95E-05|  3,605.209 
6,740.206 | 3.95E-05;  3,913.995 
7,425.351 | 3.95E-05| 4,402.936 
7,292.323 | 3.95E-05| 4,527.231 
10,679.366 | 3.95E-05|  3,434.369 
14,403.959 | 3.95E-05 | 3,750.504 
19,340.623 | 3.95E-05| 4,020.157_ 
-10,823.055 | 3.95E-05| 4,521.422 
10,959.675 | 3.95E-05| 5,078.515 
11,010.568 | 3.95E-05| 5,435.708 
10,718.151 | 3.95E-05| 5,602.052 
17,663.140 | 3.95E-05| 3,121.487 
15,522.867 | 3.95E-05| 6,085.098 _ 
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3,929.364 | 0.80055 | 32,613. 718) 33 972.965 

— 4,518.849 | 4,324.243| 4,324.243 | 10.74372 | 0.95693 | 39,313.986| 33,729.093 
— 6,090.617 | 5,893.668 | 5,893.668 | 8.84765 | 0.96766 | 35,325.579| 37,130.107 
6,483.444 | 6.876.547 | 6,483.444 | 8.04741 | 0.94283 | 36,307.287| 35,070.391 
8,054.929 | 8,054.082 | 8,054.082 | 7.49956 | 0.99989 | 37,052.672| 38,659.591 
9,037.221 | 8,839.529 | 8,839.529 | 5.59762 | 0.97812 | 28,919.107| 31 822.306 
5,107.935 | 5,107.192 | 5,107.192 | 5.24984 | 0.99985 | 19,920.947| 20,428.769 
4,322.346 | 4,322.046 | 4,322.046 | 9.49860 | 0.99993 | 32,417.597| 34,144.165 
—10,412.627 | 10,608.498) 10,412.627| 10.24927 | 0.98154 55,186.924| 59,407.586 
17,288.782 | 16,501.266! 16,501.266) 7.74731 | 0.95445 | 55,324.103| 61,054.686 
5,893.850 | 5,893.483 | 5,893.483 | 9.89697 | 0.99994 | 40,078.182| 43,611.776 
7,858.228 | 8,250.707 | 7,858.228 | 10.44928 | 0.95243 | 48,721.015| 53,629.595 
| 11,590.887 |11,589.419/ 11,589.419| 6.29947 | 0.99987 | 37,090.837| 39,404.024 
~ 2,750.663 | 2,946.921 | 2,750.663 | 9.24775 | 0.93340 | 25,031.031| 25,343.519 
-6,090.530 | 2,946.831 | 2,946.831 | 9.30052 | 0.48384 | 30,452.649| 31,825.77 
9,430.220 | 2,946.776 | 2,946.776 | 9.62474 | 0.31248 | 33,005.768| 34,477.278 
13,752.317 | 3,339.724 | 3,339.724 | 10.64554 | 0.24285 | 41,256.952| 40,744.635 
-18,271.181 | 3,339.685 | 3,339.685 | 9.24259 | 0.18278 | 34,715.244| 35,066,696 
| 2,554.207 | 6,090.189 | 2,554.207 | 7.41302 | 0.41940 | 23,498.708| 22,533.698 
2.554.194 | 9,626.315 | 2.554.194) 6.84620 | 0.26533 | 23,498.582| 23,103.157 
2,554.187 |15,126.918| 2,554.187 | 6.24797 | 0.16885 | 23,243.100| 22,690.377 
- 6,090.339 | 6,286.748 | 6,090.339 | 12.89897 | 0.96876 | 52,985.951| 57,209.402 
— 9,626.476 | 6,286.601 | 6,286.601 | 12.59431 | 0.65305 | 58,721.503| 63,494.674 
-13,359.191 | 6,286.382 | 6,286.382 | 10.54525 | 0.47057 | 52,100.845| 56,577.440 
-17,288.612 | 5,893.267 | 5,893.267 | 7.69138 | 0.34088 | 43,221.530| 40,663.544 
1,965.133 |16,701.129] 1,965.133 | 4.37760 | 0.11766 | 12,969.880} 13,360.903 
—12,968.268 | 12,377.488| 12,377.488| 7.44989 | 0.95444 | 45,388.938| 44,558.956 
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"0.49137 | 0.51185 | iis 7781 : 
82,160.612 | 47,785.441 | 0.44769 | 0.41053 | 1.07143 | 0.88148 | 0. .75626 | 0.52302 0.58161 | 12.51% | 0.056008 
87,815.649 | 46,670.951 | 0.40661 | 0.42282 | 0.93232 | 0.70567 | 0.74172 | 0.48626 | 0.53147 | 9.75% 0.039660 
86,878.151 | 45,114.351 | 0.36950 | 0.40367 | 0.90071 | 0.72487 | 0.70018 | 0.48033 | 0.51928 | 9.29% | 0.034309. 
-98,259.794 | 46,539.291 | 0.46574 | 0.39344 | 0.80260 | 0.63899 | 0.66671 | 0.44526 | 0.47364 | 8.15% | 0.037979 
—79,555.765 | 36,746.528 | 0.61548 | 0.40000 | 0.77360 | 0.60881 0.66993 | 0.43826 | 0.46190 | 8.77% 0.053963 
46,986.169 | 25,877.627 | 0.26797 | 0.43478 | 1.00538 | 0.77395 | 0.79367 | 0.50136 | 0.55075 | 12.03% 0.032228 
74,339.194 | 42,045.642 | 0.19830 | 0.45930 “4.06685 | 0 82255 | 0.86636 | 0.51618 | 0.56559 | 13.13% 0.026037 — 
163,478,247  72,420.331 | 0.20861 | 0.36340 | 0.70687 | 0.53866 0.57985 | 0.41572 | 0.44300 | 6.54% 0.019517 
184,814.183 | 68,524.345 | 0.55220 | 0.33036 | 0.55835 r 0.45079 | 0.49748 | 0.36123 | 0.37077 | 7.03% 0.038808 © 
100,189.216 | 52,125.637 | 0.30358 | 0.43529 | 0.93091 | 0.71575 | 0.77886 | 0.48212 | 0.52027 | 9.88% | 0.029982 va 
-132,018.234 | 64,942.336 | 0.43791 | 0.40623 | 0.82385 , 0.61806 | 0.68033 , 0.45658 | 0.49192 | 7.77% 0.034035 
411,258.419 | 47,113.566 | 0.74499 | 0.35417 | 0.67222 | 0.52921 | 0.56222 i 0.40200 | 0.42346 7.88% | 0.058682 
49,786.995 | 32,282.242 | 0.48538 | 0.50904 | 1.32197 | 1.02502 | 1.03782 | 0.58017 | 0.64841 , 20.32% | 0.098641 
—51,569.546 | 39,772.850 _ 0.58057 | 0.61714 | 1.51167 | 1.15742 | 1.20961 | 0.69606 | 0.77125 | 9.96% 0.057838 
§2,452.611 | 44,684.674 | 0.52290 | 0.65730 | 1.64117 | 1.21222 | 1.26627 | 0.75252 | 0.85191 | 8.21% 0.042921 
63,120.786 | 55,598.016 | 0.76730 | 0.64550 | 1.62897 | 1.20878 | 1.19377 | 0.76265 | 0.88082 | 7.80% 0.059816 — 
—§3,768.934 | 49,633.553 | 0.95714 | 0.65217 | 1.67498 | 1.17152 | 1.18338 | 0.78195 | 0.92309 | 10.96% | 0.104873 _ 
36,269.745 | 30,198.932 | 0.69697 | 0.62128 | 1.66139 | 1.29277 | 1.23968 | 0.73661 | 0.83262 | 21.87% | 0.152397 
33,715.356 | 30,457.082 | 0.74797 | 0.68524 | 1.81433 | 1.39980 | 1.37625 | 0.79169 | 0.90336 | 21.87% | 0.163547 — 
—31,416.498 | 29,794.148 | 0.17506 | 0.72224 | 1.94479 | 1.51716 | 1.48108 | 0.82920 | 0.94836 | 24.37% | 0.042671. 
132,769.396 | 69,179.766 | 0.40521 | 0.43089 | 0.91731 | 0.70258 | 0.75858 | 0.48202 | 0.52105 | 14.68% | 0.059466 _ 
130,761.309 | 75,868.277 | 0.54309 | 0.48558 | 0.99817 | 0.77257 | 0.83536 | 0.54396 | 0.58020 | 17.38% 0.094388 
108,125.774 | 67,992.707 | 0.70165 | 0.52326 | 1.06841 | 0.81868 | 0.88903 | 0.58959 | 0.62883 | 25.98% 0.182261 
74,255.167 | 51,109.352 | 0.58207 | 0.54762 | 1.17455 | 0.99327 | 0.93449 | 0.63941 | 0.68829 | 40.13% | 0.233577 
17,096.660 | 16,208.666 | 0.75862 | 0.78149 | 1.96269 | 1.57049 | 1.61784 | 0.84053 | 0.94806 | 32.09% | 0.243418 
136,152.366 | 53,773.270 | 0.33337 | 0.32727 | 0.60746 | 0.51274 | 0.50336 | 0.38117 | 0.39495 | 28.89% | 0.096314 
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0.035070 
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0.023752 
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3,086.23 
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3,075.99 
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| §,369.71 
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3,744.43 
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